Can one extract energy loss probability distributions from Raa? 
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The nuclear suppression of high transverse momentum Pt hadrons is one of the most striking 
findings in heavy-ion collision experiments. It has long been recognized that the suppression can 
be theoretically described by folding the primary parton spectrum with an energy loss probability 
distribution which is suitably averaged over the collision geometry. However, an interesting problem 
is to what degree the procedure can be inverted, i.e. given a measurement of the suppression factor 
Raa with arbitrary precision, can the probability distribution of energy loss be extracted in a 
model-independent way? In this note, we present a conceptual study of the inversion problem for 
LHC energies and demonstrate that a measurement of Raa alone is insufficient to determine the 
distribution, other observables such as 7-hadron correlations must be taken into account. 

PACS numbers: 25.75.-q,25.75.Gz 



The nuclear suppression factor Raa > defined as the 
measured spectrum of high Pt hadrons divided by the 
value seen in p-p collision scaled with the number of 
binary collisions is one of the best investigated signa- 
tures of jet quenching due to medium induced energy 
loss [l|, 0, 0, EL IE @] ■ For pions, it has been measured at 
RHIC out to 20 GeV P T by the PHENIX collaboration 
From a theory perspective, it has been realized that 
Raa needs to be computed in a probabilistic framework, 
not only because the energy loss for a given parton may 
be probabilistic (cf. e.g. but also because the initial 
vertex is not known experimentally, thus the averaging 
over the evolving geometry leads to a probability distri- 
bution of energy loss (cf. e.g. Q). 

Combining all probabilistic effects, the effect of jet en- 
ergy loss may be cast quite generally into a distribution 
(P(AE,E))t A a which describes the probability to lose 
energy AE given the initial parton energy E and an un- 
known initial vertex which is averaged over fiol . [Tlj . 

Under the assumption that the probability distribution 
for AE <C E does not strongly depend on E, one can 
reduce the computation to evaluating (P(AE))t aa - The 
medium-modified perturbative production of hadrons 
can then be computed from the (schematical) expression 
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with da-^ c r rf+x the perturbative QCD (pQCD) pro- 
duction cross section for a parton / and D^ h (z, n 2 F ) 
the fragmentation function with momentum fraction z 
at scale /sp [l2|, El (see appendix [A] and [ll| for the ex- 
plicit expression for Eq. JQ)). From this one can compute 
the nuclear modification factor Raa as 
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An intriguing question is now if it is possible to invert 
Eq. ([T]) with known Raa , assuming full knowledge of par- 
ton spectrum and fragmentation function, to determine 
(P(AE))t aa approximately in a model-independent way. 
This can be done by discretizing the integral over AE 
implicit in Eq. ([T]) which results in 
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where Raa is provided at to discrete values of Pt labelled 
Ri and (P(AE))t aa is probed at n discrete values of 
AE labelled Pj. Kij is then the calculated Raa for all 
pairs (p T ,AE J ) assuming the pQCD parton spectrum 
and fragmentation function are known (for the explicit 
expression see appendix |A"|) . 

Eq. © can be solved for the vector Pj by inversion of 
for to = n, i.e. the experimental resolution of Raa 
determines the available resolution for the energy loss 
probability density. However, in general this does not 
guarantee that the result is a probability distribution. 
Especially in the face of experimental and theoretical 
errors and finite numerical accuracy the direct matrix 
inversion may permit negative Pj which have no proba- 
bilistic interpretation. 

Thus, a more promising solution which avoids the above 
problems is to let to > n and find the vector P which min- 
imizes \\R — KP\\ 2 subject to the constraints < Pi < 1 
and J27=i Ri = 1- This guarantees that the outcome can 
be interpreted as a probability distribution and since the 
system of equations is overdetermined for to > n errors 
on individual points Ri do not have a critical influence 
on the outcome any more. 

It has been shown in [l(| that Raa for RHIC kinematics 
is very insensitive to the shape of the underlying proba- 
bility distribution. Part of the problem at RHIC is that 
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for a steeply falling parton spectrum, even a small en- 
ergy loss effectively resembles an absorption of the par- 
ton. Therefore, there is some reason to suspect that the 
inversion problem is badly conditioned. In the following, 
we study the problem at LHC kinematics where the par- 
ton spectrum is harder and larger shifts in parton energy 
can in principle be probed. 

Our setup to test the inversion is as follows: We provide 
a prediction for Raa at LHC kinematics based on [l4| 
sampled at 30 or 50 discrete points between 6 GeV and 
400 GeV momentum as input vector Ri . We solve for Pj 
at 10 or 23 discrete points between and 300 GeV energy 
loss. In addition to these discrete values, we always pro- 
vide zero (no energy loss) and 2.25 TeV (complete parton 
absorption) as possibilities and compute accordingly. 
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FIG. 1: Geometry- averaged energy loss probability density 
for LHC kinematics using a hydrodynamical model of the 
medium underlying a prediction for Raa [H[ compared with 
the probability density extracted from the predicted Raa via 
numerical inversion for different numbers m,n (see text). 

The resulting solution to the vector Pj interpreted as a 
probability density is shown in Fig. [1] and compared with 
the (computed) input energy loss probability distribution 
underlying the prediction of Raa in UM- While details 
of the distribution depend on the values of m and n, 
the main features of the extracted distribution are rather 
stable: There is some strength close to zero energy loss, 
the distribution peaks between 50 and 70 GeV energy 
loss and the strength in the absorption channel (not seen 
in Fig.[T]) is about 0.35, i.e. roughly a third of partons is 
absorbed. However, the most striking observation is that 
the distribution found from the numerical inversion does 
not resemble the input distribution closely. 
One may ask if the extracted probability distribution is 
able to describe Raa when inserted back into Eq. (JTJ) . 
This is shown for m = 30, n = 25 in Fig. [5] where we 
compare the original prediction of Raa with the one re- 
covered from the extracted probability density. As can 
easily be verified, the results cannot be distinguished, 
i.e. two very different probability distributions yield an 
indistinguishable Raa ■ 

We are thus faced with a rather asymmetric problem — 
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FIG. 2: Raa as predicted for LHC kinematics using a hy- 
drodynamical model for the medium evolution [lj] compared 
with Raa computed from the probability distribution ex- 
tracted by numerical inversion of Raa (see text). The two 
curves are virtually identical, hence a shifted recovered curve 
is provided for comparison. 



the fact that Raa is rather insensitive to the shape of the 
probability distribution prevents the extraction of any- 
thing but gross features of the distribution given only a 
measurement of Raa- Since in reality neither the pQCD 
spectrum of partons nor the fragmentation function is 
free of uncertainties and since a measured Raa will have 
experimental errors, the strategy to invert Raa to find 
the energy loss probability distribution is simply not fea- 
sible in practice. 
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FIG. 3: Extracted energy loss probability density from the 
prediction of Raa at LHC under the assumption that the 
EKS98 nuclear parton distribution or the NPDF nuclear par- 
ton distribution is used to compute the folding kernel Kij (the 
original prediction is based on EKS98). 

We illustrate the the uncertainty associated with the 
knowledge of the pQCD kernel K%j in Fig. [3J Here, we 
solve Eq. ([3]) using two different sets of nuclear parton 
distributions to compute Ky, i.e. the EKS98 set [l5| 
and the NPDF set [l6| (note that the prediction for Raa 
used as input is based on EKS98). It is clearly seen that 
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the (rather small) uncertainty in the choice of the par- 
ton distribution function maps into visible changes in the 
extracted probability distribution. While gross features 
and typical energy loss scales are still similar, the differ- 
ences in detail are substantial. 

Thus, while experimental data on Raa are valuable in 
constraining energy loss models, due to the badly con- 
ditioned inversion problem Raa per se does not even 
measure an energy loss probability distribution, much 
less a density-distribution of matter (which can only in 
a model-dependent way be linked to the energy loss dis- 
tribution) and even less a single parameter character- 
izing the density. Turning to other, more differential 
observables (such as 7-hadron correlation measurements 
in which the energy of a quark is known before energy 
loss |l0() is necessary to achieve the first step towards a 
model-independent extraction of medium properties from 
hard probes. 



APPENDIX A: CALCULATING K i3 

In this appendix, we describe how to compute the inver- 
sion kernel as a function of hadron momenta P T and 
energy loss AP 7 from perturbative QCD (pQCD). 

Let A and B stand for the colliding objects (protons or 
nuclei) and be yi(2) the rapidity of parton k(l). The 
distribution function of a parton type i in A at a mo- 
mentum fraction X\ and a factorization scale Q ~ px is 
fi/A{x\,Q 2 )- The distribution functions are different for 
the free protons [H, 0j[ and protons in nuclei . 
The fractional momenta of the colliding partons i, j are 
given by xi <2 = ^ (exp[±j/ x ] + exp[±y 2 ]). 

Inclusive production of a parton flavour / at rapidity y f 
is found by integrating over either y\ or yi and summing 
over appropriate combinations of partons, 
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where the summation (ij),(kl) runs over pairs 
99i 9Q.i 59; 111 <Z9i QQ an d 9 stands for any of the quark 
flavours u,d,s. Expressions for the pQCD subprocesses 
da — (s,i,ii) as a function of the parton Mandelstam 

variables s,i and u can be found e.g. in [l7j |. 
For the production of a hadron h with mass M, trans- 
verse momentum Pt at rapidity y and transverse mass 
niT = J M 2 + P^ from the parton /, let us introduce 
the fraction z of the parton energy carried by the hadron 
after fragmentation with z = Eh/ Ef. Assuming collinear 



fragmentation, the hadronic variables can be written in 
terms of the partonic ones as 



niT cosh y — zpx cosh y f and tot sinh y = Pt sinh y f . 

(A2) 

Thus, the hadronic momentum spectrum arises from the 
partonic one by folding with the probability distribution 
Df^h(z, fxi) to fragment with a fraction z at a scale (if ~ 
Pt as 



da AB^h+X r da AB ^f+ X 



r d(J AB^j+X rl 

dP^dy ~ E J d Pr d Vf dp2 dyf J dzD f ^ h (z,fi 2 f )6(m T ~ M T (p T ,y f ,z))S(y-Y(p T ,y f ,z)) (A3) 
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with and 

Y(p T , y f ,z) = arsinh f^- sinht//^ . (A5) 

The lower cutoff z m i n = 2 ^= L cosh y arises from the fact 
that there is a kinematical limit on the parton momen- 
M T (pT,yf,z) = (zpt) 2 + M 2 tanh 2 y/ (A4) turn; it cannot exceed \fs/{2 cosh yf) and thus for given 
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hadron momentum there is a minimal z corresponding to 
fragmentation of a parton with maximal momentum. 

I 



As a side remark, in order to obtain the modified spec- 
trum for a probabilistic energy loss instead of a fixed 
value of AE, Eq. (|A6|) has to be integrated over the 
probability of energy loss J cIAE(P(AE)}t aa - Insert- 
ing Eq. into Eq. (KM instead of dg ^~y + * yields 

i 
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